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AX  INVESTIGATION  OF  THE  IGNITION  DELAY  TIMES 
FOR  PROPYLENE  OXIDE-OXYGEN- NITROGEN  MIXTURES 


ABSTRACT 


A reflected  shock  tube  technique  has  been  used  to  measure  the  ignition 
delay  time  for  propylene  oxide-oxygen-nitrogen  mixtures  in  the  temperature 
range  of  900  to  1250°K.  Mixture  equivalence  ratios  ranging  from  <p  = 0.8  to 
$ = 1.25  were  used.  Ignition  delay  times  were  measured  by  using  both  streak 
schlieren  and  pressure  gauges  mounted  on  the  back  wall  of  the  shock  tube.  A 
regression  analysis  indicated  that  the  delay  time  multiplied  by  the  pro- 
pylene oxide  concentration  to  the  .8  power  and  the  oxygen  concentration  to 
the  -1.2  power  yielded  a best  fit  when  plotted  as  a logarithm  against  the 
reciprical  temperature.  This  regression  analysis  also  indicated  that  the 
nitrogen  concentration  had  no  effect  on  the  ignition  delay  time.  The  ap- 
parent activation  energy  for  the  reaction  was  found  to  be  approximately 
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I.  INTRODUCTION 


In  the  past  few  years  an  increasing  number  of  experimental  studies 
have  been  conducted  on  modeling  the  ignition  delay  times  involved  in  fuel- 
air  explosions.  Reasons  for  these  studies  can  be  attributed  to  the  search 
for  alternative  fuel  sources,  to  the  prevention  of  explosions  occurring 
after  accidents  involving  highly  combustible  fuels,  and  to  the  development 
of  fuel-air  explosion  weaponry. 

A literature  search  revealed  that  recent  research  in  fuel-air  ex- 
plosions has  been  conducted  using  propylene  oxide  [1,2].  Yet  there  appears 
to  be  a lack  in  the  literature  of  any  ignition  delay  models  involving  pro- 
pylene oxide  which  might  aid  in  analyzing  these  results.  An  ignition  study 
was  conducted  by  Hikita  and  Yoneda  using  a shock  tube  [3].  An  abstract  of 
their  report  [4]  shows  that  in  propylene  oxide-oxygon-argon  mixtures  with 
fuel  equivalence  ratios  between  0.353  and  2.67,  the  rate  controlling  re- 
action's activation  energy  is  1.76  x 105  joule/mole  (42.1  Kcal/mole).  No 
mention  was  made  of  the  concentration  dependency  of  the  overall  rate. 

In  this  investigation  the  activation  energy  and  the  concentration 
dependencies  of  the  rate  controlling  reaction  in  the  ignition  of  propylene 
oxide-oxygen-nitrogen  mixtures  are  determined  by  correlating  the  experi- 
mental data  to  the  Arrhenius  relation:  t^  = Aexp(E/RT)  [C_H^0]a [021^^2]° . 
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II.  EQUIPMENT 

A.  University  of  Illinois  Shock  Tube 

The  ignition  study,  described  in  this  investigation,  was  conducted 
using  the  buffer  chamber  and  the  low  pressure  chamber  of  the  Aeronautical 
and  Astronautical  Engineering  Department's  rectangular  shock  tube.  A brief 
description  of  the  shock  tube  follows;  a more  detailed  description  can  be 
found  in  Reference  5. 

The  buffer  chamber,  used  as  the  high  pressure  (compression)  chamber, 
is  4.04  meters  long  and  has  an  internal  cross  section  of  0.1016  x 0.381  meters. 

The  low  pressure  chamber  is  separated  into  two  sections:  an  expansion 
section  8.64  meters  long  and  a test  section  0.864  meters  long.  Both  sections 
have  an  internal  cross  section  of  0.1016  x 0.381  meters.  The  test  section 
contains  two  0.14  meter  circular  windows  for  optical  investigation  of  the 
test  section's  flow  field.  Mounted  inside  the  test  section,  overlapping  the 
circular  windows,  a 50  mm  thick  aluminum  block  forms  the  rear  wall  upon  which 
the  incident  shock  reflects. 

The  test  section  and  the  buffer  section  are  mounted  on  rollers  and  are 
connected  to  the  expansion  section  by  means  of  hydraulic  cylinders  and  con- 
necting rods. 

B.  Pressure  Transducer  System 

Previous  detonation  studies  conducted  with  the  Aeronautical  and  Astro- 
nautical  Engineering  Department's  rectangular  shock  tube  [5-7]  have  used  a 
combination  optical  interferometer  and  light  field  schlieren  system  to 
measure  the  incident  shock  velocity,  the  reflected  shock  velocity,  the  den- 
sity ratios  across  the  shocks,  and  the  ignition  delay  time.  In  this  study 
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and  its  immediate  predecessor  [8],  the  incident  shock  velocity  and  the  igni- 
tion delay  time  were  measured  using  a pressure  transducer  system  composed 
of  two  Kistler  103A  pizoelectric  transducers,  a Kistler  503  charge  ampli- 
fier, a Kistler  504  dual  mode  amplifier,  a Tektronix  555  dual  beam  oscillo- 
scope and  camera,  and  a Potter  1.6  megacycle  counter  chronograph.  A block 
diagram  of  the  system  is  shown  in  Fig.  1. 

One  transducer  is  mounted  in  the  shock  tube's  side  wall  (side  wall 
transducer)  0.624  meters  upstream  of  the  aluminum  block.  As  the  incident 
shock  passes  the  side  wall  transducer,  the  transducer  produces  a signal  which 
triggers  the  oscilloscope's  lower  beam  and  starts  the  event  counter.  The 
second  transducer  is  mounted  in  the  center  of  the  aluminum  block  (rear  wall 
transducer);  the  signals  produced  from  this  transducer  are  recorded  on  the 
oscilloscope's  lower  and  upper  beams.  When  the  incident  shock  strikes  the 
aluminum  block,  the  rear  wall  transducer  produces  a signal  which  triggers 
the  oscilloscope's  upper  beam  and  stops  the  event  counter.  The  oscilloscope's 
traces  are  recorded  on  3090  ASA  Polaroid  film. 

Example  pictures  are  shown  in  Figs.  2 and  3.  The  oscilloscope's  upper 
beam  traces  were  set  at  sweep  speeds  varying  from  10  ysec/cm  to  200  ysec/cm; 
this  allowed  reducible  pressure-time  histories  of  the  ignition  processes  to 
be  recorded.  The  oscilloscope's  lower  beam  traces  were  set  at  a sweep  speed 
of  200  ysec/cm;  this  allowed  pressure-time  histories  from  the  incident  shocks 
passing  the  side  wall  transducer  through  the  detonation  to  be  recorded. 


Mixture  Gases  and  Gauges 


The  propylene  oxide  used  in  this  investigation  was  produced  by  East- 


man, Rochester,  Wv.  The  gases  used  were  high  purity  helium,  extra  dry  oxygen 
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and  prepurified  nitrogen  bottled  by  Linde.  Two  Wallace  and  Tiernan  (800  and 
1600  Torr)  precision  dial  manometers  were  used  to  measure  the  pressures  used 
in  mixing,  and  to  measure  the  pressures  of  the  low  pressure  chamber  and  the 
high  pressure  chamber  during  an  experiment. 
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III.  PROCEDURE 

A.  Mixing 

Six  hours  preceding  a sequence  of  tests,  the  mixing  chamber  and  the 
bleedoff  chamber  were  evacuated  below  a pressure  of  1.33  pascals  and  isolated. 
Using  a syringe,  liquid  propylene  oxide  was  injected  into  the  evacuated  mix- 
ing chamber  through  a compressed  neoprene  diaphragm.  Using  partial  pressures, 
the  excess  propylene  oxide  was  bled  off;  the  appropriate  amounts  of  nitrogen 
and  oxygen  were  then  added  to  obtain  the  desired  mixtures  listed  in  Table  1. 
For  safety  purposes,  the  total  pressure  of  the  test  mixture  in  the  mixing 
chamber  was  maintained  below  5.45  x 104  pascals  (410  Torr) . 


B.  Test  Run 

Preparation  for  a test  run  began  with  the  insertion  of  a diaphragm  be- 
tween the  high  pressure  (buffer)  chamber  and  the  low  pressure  chamber  of  the 
shock  tube;  76.2  pm  acetate  diaphragms  were  used  for  driving  pressures  under 
1.15  x 10s  pascals  while  127  pm  acetate  diaphragms  were  used  for  driving 
pressures  over  1.15  x 105  pasclas.  The  shock  tube  was  then  evacuated;  during 
this  time  the  electronics  were  switched  on  and  allowed  to  warm  up.  When  the 
shock  tube  pressure,  measured  by  an  Alphatron  Gauge,  reached  a pressure  lower 
than  10  pasclas,  the  high  pressure  chamber  was  isolated  and  filled  with  helium. 
As  the  high  pressure  chamber's  pressure  neared  three  quarters  the  driving  pres- 
sure, the  low  pressure  chamber  was  isolated  from  the  vacuum  pump  and  filled 
with  the  test  mixture  to  a pressure  of  1.33  x 103  pascals  (10  Torr) . After 
allowing  enough  time  for  the  test  mixture  to  attain  equilibrium,  the  tempera- 
ture of  the  text  mixture  was  recorded  using  a thermometer  placed  in  the  top  of 
the  shock  tube  near  the  test  section.  Just  before  rupturing  the  diaphragm, 
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all  equipment  was  rechecked  and  set;  the  oscilloscope's  camera  shutter  was 
then  opened. 

After  the  shot  was  taken  the  oscilloscope's  camera  shutter  was  closed, 
the  camera  film  developed,  the  Potter  event  counter  time  was  recorded  and 
the  shock  tube  was  cleaned  and  prepared  for  another  experiment. 

C.  Data  Reduction 

The  letter  event  counter  time  and  the  test  mixture's  initial  tempera- 
ture, pressure  and  composition  were  used  as  input  for  the  computer  program 
SH3.F4.  The  program  initially  calculates  the  incident  shock  velocity  by 
dividing  the  distance  between  the  side  wall  and  the  rear  wall  transducers  by 
the  time  measured  on  the  Potter  counter.  Using  the  incident  shock  velocity 
and  the  other  inputted  data,  the  program  then  procedcs  to  calculate  and  out- 
put the  temperature  and  the  concentrations  of  propylene  oxide,  oxygen  and 
nitrogen  behind  the  reflected  shock.  A description  of  the  program  is  given 
in  Appendix  A. 

There  were  two  nodes  of  ignition  observed  in  this  investigation,  an 
occurrence  which  has  been  studied  in  H^-O^  detonation  by  Meyer  and  Oppenheim 
[9].  Figure  2 represents  a strong  ignition;  on  the  upper  beam  trace  the  in- 
stantaneous appearance  of  the  strong  pressure  spike  is  associated  with  a 
one-dimensional  detonation  initiation.  Figure  3 represents  a mild  ignition; 
on  the  upper  beam  trace  the  appearance  of  the  mild  pressure  spike  is  associ- 
ated with  the  formation  of  distinct  flame  kernels,  the  appearance  of  a rela- 
tively constant  pressure  following  the  mild  pressure  spike  is  associated  with 
the  growtli  of  the  flame  kernels,  and  the  appearance  of  the  second  pressure 
spike  is  associated  with  the  formation  of  a detonation  wave.  The  ignition 
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delay  time  was  calculated  from  the  distance  between  the  triggering  and  the 
first  pressure  spike  recorded  on  the  upper  beam  trace,  the  first  pressure 
spike  representing  the  beginning  of  ignition  in  both  strong  and  mild  igni- 
tions. 

Presented  in  Tables  2-5  are  the  reflected  shock  properties  calculated 
from  SH3.F4  and  the  ignition  delay  times  measured  from  the  oscilloscope's 
upper  beam  traces. 

Plots  of  the  ignition  delay  times  versus  the  reciprocal  of  the  re- 
flected shock  temperatures  for  each  mixture  are  presented  in  Figs.  4-7. 
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TABLE  2 

KINETIC  DATA  FOR  SHOCK  TUBE  RUNS 

Fuel  Equivalence  Ratio  =1.0 
C,H^0  Percentage  = 18% 


103/T 

°K"  1 

(C3H6°) 

10*2  K-mol/m3 

C°2) 

10-2  K-mol/m3 

(n2) 

10~2  K-mol/m3 

Initiation 
Delay  Time 

10" 6 sec 

0.955 

0.23SS1 

0.95523 

0.13267 

230.0 

0.874 

0.26965 

1.07860 

0.14981 

44.0 

0.85S 

0.27613 

1 . 10450 

0.15340 

26.5 

0.928 

0.24852 

0.99408 

0. 13807 

110.0 

1.035 

0.21102 

0.84406 

0.11723 

1180.0 

0.888 

0.27237 

1.08950 

0.15132 

65.0 

0.943 

0.25127 

1.00510 

0.13959 

152.5 

0.968 

0.24210 

0.96841 

0. 1345 

460.0 

0.982 

0.23651 

0.94604 

0. 13139 

520.0 

0.980 

0.25964 

0.95856 

0.13313 

435.0 

0.993 

0.23283 

0.93130 

0.12955 

675.0 

0.867 

0.27617 

1.10470 

0.15343 

48.0 

0.885 

0.26665 

1.06650 

0.14813 

58.0 

1.035 

0.21612 

0.S6450 

0. 12007 

1140.0 

1.024 

0.2230 

0.88921 

0.12550 

1060.0 

0.966 

0.23862 

0.95447 

0. 13256 

205.0 

0.920 

0.25316 

1.01260 

0. 14064 

88.0 

0.902 

0.259S1 

1.03920 

0.14434 

64.0 

0.932 

0. 24"29 

0.98914 

0.13738 

110.0 

0.973 

0.23204 

0.92814 

0.12891 

227.5 

0.996 

0. 23665 

0.94659 

0. 13147 

350.0 

0.965 

0.23596 

0.94384 

0.13109 

192.5 

0.922 

0. 2569" 

1 .02790 

0.14276 

162.5 

0.970 

0.23329 

0.93314 

0.12960 

255.0 

0.920 

0. 24937 

0.99750 

0.13854 

98.0 

1.037 

0. 20~84 

0.83139 

0.11547 

1050.0 

1.026 

0.21003 

0.84014 

0.11669 

800.0 

0.969 

0. 22543 

0.90174 

0.12524 

250.0 

0.904 

0. 24‘12 

0.98850 

0.13729 

77.0 

0.992 

0. 22553 

0.90133 

0.12518 

420.0* 

0.995 

0.22901 

0.91605 

0.12723 

435.0* 

0.975 

0.23"97 

0.94988 

0.13193 

300.0* 

1.022 

0. 20940 

0.83758 

0.11633 

490.0* 

Mild  ignitions. 


KINETIC  DATA  FOR  SHOCK  TUBE  RUNS 


Fuel  Equivalence  Ratio  = 1.25 
C_HgO  Percentage  = 21.429% 


io3/t 

°K~ 1 

(C3H6°) 

10  2 K-mol/n3 

C°2) 

10  2 K-mol/m3 

(n2) 

10  2 K-mol/m3 

Initiation 
Delay  Time 

10  6 sec 

1.035 

0.28592 

0.90S54 

0.13250 

880.0 

0.966 

0.31635 

1.01230 

0.14763 

195.0 

0.936 

0.33678 

1.07770 

0.15717 

116.0 

0.985 

0.31091 

0.99493 

0.14509 

235.0 

1.008 

0.29778 

0.95289 

0.13896 

355.0 

1.064 

0.27134 

0.86828 

0.12662 

1070.0 

0.923 

0.34419 

1.10140 

0.16062 

72.0 

0.955 

0.32720 

1.04700 

0.15269 

148.0 

1.061 

0.27730 

0.88737 

0.12941 

1000.0 

0.993 

0.30550 

0.97760 

0.14257 

280.0* 

1.065 

0.28915 

0.92528 

0.13494 

1040.0 

1. 013 

0.31027 

0.99288 

0. 14479 

310.0 

0.954 

0. 33586 

1 .07480 

0.15674 

134.0 

1.019 

0.29709 

0.95068 

0.13864 

360.0 

1.027 

0.29146 

0.93266 

0. 13601 

555.0 

0.984 

0.31600 

1.01120 

0.14747 

245.0 

0.955 

0.32086 

1.02670 

0.14973 

160.0 

1.023 

0.28470 

0.90934 

0.13261 

590.0 

1 . 033 

0.28817 

0.92214 

0.13448 

875.0 

0.964 

0.30695 

0.98224 

0.14324 

195.0 

1.041 

0.27162 

0.86919 

0.12676 

960.0 

1.014 

0.29147 

0.93269 

0.13602 

340.0 

0.983 

0.30566 

0.97810 

0.14264 

240.0 

0.933 

0.52S91 

1.05250 

0.15349 

117.0 

1.038 

0.27724 

0.88716 

0.12938 

970.0 

1.018 

0.30183 

0.97584 

0.14085 

350.0 

0.S94 

0.35912 

1 . 14920 

0.16759 

93.0 

0.900 

0.24234 

1.09550 

0.15976 

99.0 

0.895 

0.35469 

1.13500 

0. 16552 

94.0 

Mild  ignitions 


14 


TABLE  4 

KINETIC  DATA  FOR  SHOCK  TUBE  RUNS 

Fuel  Equivalence  Ratio  = .8 
C^H^O  Percentage  =15% 


1 o3/t 

°K-J 

CC3H6°) 

10  2 K-mol/m3 

C°2) 

10~2  K-mol/m3 

(n2) 

10  2 K-mol/m3 

Initiation 
Delay  Time 

10"2  sec3 

1.040 

0.15971 

0.79854 

0.10647 

1210.0 

0.885 

0.19915 

0.99576 

0.13277 

62.0 

0.91S 

0.1S882 

0.94410 

0.12588 

114.0 

1.013 

0.16334 

0.81669 

0. 10889 

635.0 

0.985 

0.16536 

0.82678 

0.11024 

355.0 

0.976 

0. 19431 

0.97153 

0.12954 

48.0 

0.935 

0.17727 

0.88637 

0.11818 

155.0 

1.011 

0.15711 

0.78555 

0.10474 

810.0 

0.892 

0.18983 

0.94814 

0.12655 

59.0 

1.050 

0.15410 

0.77049 

0.10273 

1080.0 

0.890 

0.18836 

0.94182 

0.12558 

61.0 

0.950 

0.17263 

0.86317 

0.11509 

190.0 

0.906 

0.18753 

0.93766 

0.12502 

75.0 

0.922 

0. 18243 

0.91215 

0.12162 

125.0 

0.962 

0.17077 

0.85386 

0.11385 

240.0 

0.884 

0.19545 

0.97727 

0.13030 

44.0 

0.929 

0.18727 

0.93636 

0.12485 

136.0 

0.932 

0.1896S 

0.94842 

0.12646 

104.0 

0.954 

0.18248 

0.91239 

0.12165 

195.0 

0.968 

0.17722 

0.8S610 

0.11815 

260.0 

0.991 

0.17171 

0.85856 

0.11447 

435.0* 

0.999 

0.16085 

0.80474 

0.10730 

480.0* 

1.006 

0.15936 

0.79681 

0. 10624 

490.0* 

1.020 

0. 15575 

0.77863 

0. 10382 

690.0* 

Mild  ignitions 
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TABLE  5 

KINETIC  DATA  FOR  SHOCK  TUBE  RUNS 

Fuel  Equivalence  Ratio  =1.0 
C,H^O  Percentage  =16% 


io3/t 

°k”  1 

CC3H6°) 

10'2  K-mol/m3 

(°2) 

10  2 K-mol/m3 

(n2) 

10’2  K-mol/m3 

Initiation 
Delay  Time 

10  6 sec 

0.903 

0.21249 

0.84997 

0.26562 

66.0 

0.902 

0.21322 

0.S5288 

0.26653 

62.0 

0.943 

0.19970 

0.79879 

0.24962 

175.0 

0.915 

0.21342 

0.85366 

0.26679 

83.0 

0.951 

0.20127 

0.80507 

0.25159 

172.5 

0.974 

0.19288 

0.77153 

0.24110 

250.0 

1.032 

0.17558 

0.70234 

0.21948 

1020.0 

0.991 

0.18798 

0.75191 

0.23497 

332.5 

0.926 

0.20448 

0.81792 

0.25560 

122.0 

0.901 

0.21111 

0.84445 

0.26389 

75.0 

0.982 

0.18731 

0.74922 

0.23413 

277.5 

0.953 

0.19262 

0.77048 

0.24077 

150.0 

0.898 

0.21266 

0.85063 

0.26582 

72.0 

0.925 

0.20538 

0.81350 

0. 25422 

132.0 

0.978 

0.18711 

0.74845 

0.23389 

230.0 

1.014 

0.17787 

0.71 149 

0.22234 

425.0* 

1.001 

0.18421 

0.736S3 

0.23026 

320.0* 

1.023 

0.17825 

0.71300 

0.22281 

500.0* 

1.017 

0. 17654 

0.70617 

0. 22068 

680.0* 

1.029 

0. 17115 

0.68462 

0.21394 

730.0* 

0.987 

0.18604 

0.74416 

0.23255 

300.0* 

1.024 

0.17150 

0.68599 

0.21437 

595.0 

1.029 

0. 16981 

0.67924 

0.21226 

885.0 

0.946 

0.1908S 

0.76352 

0.19088 

125.0 

0.949 

0. 18787 

0.75149 

0.23484 

142.5 

Mild  ignition 
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1000/T  (°K  ’) 


IGNITION  DELAY  TIME  VS.  10 3 /REFLECTED  SHOCK  TEMPERATURE 


FIGURE  TV 
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1000/T  (°K  *) 


IGNITION  DELAY  TIME  VS.  IOVrEFLECTF.D  SHOCK  TEMPERATURE 


FIGURE  VI 


1000/T  (°K  *) 


IGNITION  DELAY  TIME  VS.  103/REFLECTED  SHOCK  TEMPERATURE 


FIGURE  VII 
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IV.  RESULTS  AND  DISCUSSION 


The  goal  of  this  investigation  was  to  develop  a model  relating  the 
ignition  delay  time  to  the  temperature  and  the  gas  coi.  jntrations  in  pro- 
pylene oxide-oxygen-nitrogen  detonations.  Using  the  multiple  regression 
routine  in  Appendix  B,  the  Arrhenius  rate  relation  for  calculating  the  ig- 
nition delay  time  in  seconds  was  found  to  be: 


t.  = 6.203  x 10 


- i 2 


exp< 


flS.924  x 103 
1 T 


1C3h6°1 


'0.812 


[02]  1 - 22  5 [n2] 


0.003 


The  activation  energy,  using  a 90°^  confidence  interval,  was  found  to  be 
(1.573  ±0.209)  x 105  Joules/mole. 

Two  checks  were  made  to  determine  the  reliability  of  these  results. 
First,  theory  indicates  that  nitrogen,  because  of  its  high  dissociation 
energy,  should  act  as  an  inert  in  the  temperature  ranges  of  this  detonation 
study.  As  an  inert,  nitrogen  should  have  a concentration  dependency  expo- 
nent near  zero;  the  experimental  concentration  dependency  exponent  0.003  is 
very  satisfactory.  Second,  the  four  mixtures'  kinetic  data  were  reduced  to 
a single  line  by  plotting  t^ [C.H^O] 0 * 8 1 2 [0^]  1 ' 2 2 5 [N9] ~° ‘ 0 0 3 versus  103/T 
in  F’ig.  8.  The  fit  along  the  least  squares  was  satisfactory. 

The  experimental  Arrhenius  rate  relation  indicates  that  propylene 
oxide  concentration  increases  promote  detonation  while  oxygen  concentration 
increases  inhibit  detonation.  In  direct  contrast  Cj-C^.  alkane  research  con- 
ducted by  Burcat,  Scheller,  and  Lifshitz  [10]  indicates  that  alkane  concen- 
tration increases  inhibit  detonations. 

Approximately  4 0°s  of  the  ignition  within  the  temperature  range  970  - 
1025°K  were  mild  ignitions.  In  an  attempt  to  explain  the  differences 
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between  the  mild  ignitions  and  the  strong  ignitions,  reference  is  made  to 
the  ignition  studies  by  Meyer  and  Oppenheim  [9],  Voevodshy  and 

Solouklin  [11],  and  dohen  and  Larsen  [12].  In  these  studies  mild  ignitions 
were  characterized  by  the  appearance  of  distinct  flame  kernels  during  the 
initial  stages  of  ignition.  The  mild  ignitions  also  occurred  near  the  second 
explosion  limit,  where  chain  reaction  sensitivity  to  gas-phase  recombinations 
increases.  In  this  investigation  mild  ignitions  are  assumed  to  occur  near 
the  second  explosion  limit.  Because  of  the  chain  reactions  sensitivity  to 
gas-phase  recombination,  distinct  flame  kernels  appear  instead  of  a strong 
ignition.  The  rate  of  heat  released  from  these  flame  kernels  are  too  small 
to  generate  a strong  pressure  pulse,  thus  the  weak  pressure  pulse  found  on 
the  upper  beam  trace  in  Fig.  3. 
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APPENDIX  A 

The  program  SH3.F4  is  a revised  version  of  the  Fortran  Program  SH2.F4 
[8],  modified  to  handle  propylene  oxide-oxygen-nitrogen  mixtures.  The  pro- 
gram inputs  are:  the  test  mixture-test  run  number,  the  test  mixture's 
partial  pressures,  the  shock  tube  mixture's  initial  temperature  and  initial 
pressure,  and  the  Potter  counter  time  measurements.  The  SH3.F4  program 
outputs  are:  the  fuel  equivalence  ratio,  the  fuel  percentage,  the  average 
molecular  weight,  the  specific  heat  ratio,  the  incident  Mach  number  and 
shock  conditions,  and  the  kinetic  parameters. 

The  equations  used  in  solving  the  incident  and  reflected  shock  condi- 
tion assume  a one-dimensional,  inviscid,  steady  state,  non-reactive,  ideal 
gas  flow  system.  The  shock  equations  used  in  SH2.F4  and  SH3.F4  programs 
are  derived  in  Reference  8.  Essentially  the  program  is  an  iteration  scheme 
which  assumes  a temperature  behind  a shock  wave,  calculates  a shock  velocity, 
and  compares  this  value  with  known  conditions.  For  convergence  the  cal- 
culated square  of  the  incident  shock  velocity  must  equal  the  square  of  the 
measured  incident  shock  velocity  within  one  part  in  1 04 , and  the  calculated 
velocity  at  the  back  wall  must  be  less  than  one  part  in  104  of  the  incident 
flow  velocity. 
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APPENDIX  B 


The  multiple  regression  routine  in  a Soupac  Program  run  on  the  Uni- 
versity of  Illinois  360  computer  system.  The  program  is  based  on  the  matri- 
cized  linear  regression  model  presented  in  Neter  and  IVasserman  [13]. 

A general  linear  regression  model  can  be  expressed  as: 


y.  = a + 
J x o 


ax. 

1 1,1 


+ a_x.  . 

2 1,2 


.a  , X . , + 

p-1  i,p-l 


e. 

l 


e^  is  the  random  error  terra.  Rewritten  in  matrix  form,  the  model  becomes: 

Y = X A + E (1) 

nxl  nxp  pxl  nxl 


The  least  squares  estimators  of  A are: 

A = (X'X)"1X'Y 

The  sample  variance  is  equal  to: 

Sz  = (Y'Y  - AX' Y)/n-p 


(2) 


(3) 


In  this  ignition  delay  time  study,  the  multiple  regression  model  chosen 
was  the  Arrhenius  rate  relation: 


*ij  ■ A »p[Eo/RT.nC3H(0]“  [02lJ  [N2]'  c. 

tj  = ignition  delay  time  (sec) 

A = pre-exponential  constant 

Eq  = activation  energy  (Joules/mole) 

T = reflected  shock  temperature  (°K) 

R = universal  gas  constant  (8.314  Joules/mole-°K) 
a,b,c  = exponents 


(4) 


— 
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[ ] = concentrations  (K-moles/meter 3) 
e = random  error  term 


Linearization  of  the  regression  model  yields: 

Ant„  = AnA  + (Eq/R)  (1000/T. ) + a AntC.^0]  + b + c An[N2] . 

(5) 

Essentially  the  Soupac  Program  solves  the  matrix  equations  2 and  3 
after  equation  5 is  put  in  the  form  of  equation  1.  The  Soupac  Program  cal- 
culations for  the  least  squares  estimation  of  the  parameters  in  equation  5 
are: 


-25.806 

E /R 

0 

18.925 

a 

= 

- 0.812 

6 

1 . 225 

c 

0.003 

For  n = 111  and  p = 5,  the  sample  variance  (s2)  equaled  0.05. 

A 90%  confidence  interval  was  calculated  for  Eo/R  using  the  F dis- 
tribution table  in  Reference  13  and  the  confidence  interval  relation: 


|t/R  ■ eo/r|  < cf.90>  ,,  106),/J  S 


(O/s) 2 = [0  1 0 0 0]  (X’X)'1 
'll 


The  confidence  interval  for  E /R  is  18.925  ±2.52  °K 
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